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Abstract 

This Letter reports a measurement of the CP violation observables Sj^ K o and 
Cj^ K o in the decay channel B° — > J/tpK® performed with l.Ofb -1 of pp collisions 
at -y/s = 7TeV collected by the LHCb experiment. The fit to the data yields 
Sj /i>K o = 0.73 ± 0.07 (stat)± 0.04 (syst) and C J/%pK o = 0.03±0.09 (stat) ±0.01 (syst). 
Both values are consistent with the current world averages and within expectations 
from the Standard Model. 
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1 Introduction 



The source of CP violation in the electroweak sector of the Standard Model (SM) is the 
single irreducible complex phase of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing 
matrix [lj[2j. The decay B° — > J/ipK® is one of the theoretically cleanest modes for the 
study of CP violation in the B° meson system. Here, the B° and B° mesons decay to a 
common CP-odd eigenstate allowing for interference through B°-B° mixing. 

In the B° system the decay width difference AIT^ between the heavy and light mass 
eigenstates is negligible. Therefore, the time dependent decay rate asymmetry can be 
written as |3j|4) 



Here B°(t) and B°(t) are the states into which particles produced at t = as B° and 
B° respectively have evolved, when decaying at time t. The parameter Am^ is the mass 
difference between the two B° mass eigenstates. The sine term results from the interference 
between direct decay and decay after B°-B° mixing. The cosine term arises either from 
the interference between decay amplitudes with different weak and strong phases (direct 
CP violation) or from CP violation in B°-B° mixing. 

In the SM, CP violation in mixing and direct CP violation are both negligible in 
B° — > J/ipK® decays, hence Cj/^ K o ~ 0, while Sj^ K o rs sin 2/3, where the CKM angle 
(3 can be expressed in terms of the CKM matrix elements as arg |— V cd V* b /V td V t * b \. It 
can also be measured in other B° decays to final states including charmonium such as 
J/ipK®, J/ipK*°, ip(2S)K(*>°, which have been used in measurements by the BaBar and 
Belle collaborations |5|[6). Currently, the world averages are Sj/^ K o = 0.679 ± 0.020 and 
CjAMS* = 0-005 iO.oVfg. 

The time-dependent measurement of the CP parameters Sj^ K o and Cj^ K o requires 
flavour tagging, i. e. the knowledge whether the decaying particle was produced as a B° 
or a B° meson. If a fraction u of candidates is tagged incorrectly, the accessible time- 
dependent asymmetry Aj/^ K o(t) is diluted by a factor (1 — 2ui). Hence, a measurement of 
the CP parameters requires precise knowledge of the wrong tag fraction. Additionally, the 
asymmetry between the production rates of -B and B° has to be determined as it affects 
the observed asymmetries. 

In this Letter, the most precise measurement of Sj/^ K o and Cj^ K o to date at a hadron 
collider is presented using approximately 8200 flavour-tagged B° — > J/ip K® decays. 

2 Data samples and selection requirements 

The data sample consists of 1.0 ftT 1 of pp collisions recorded in 2011 at a centre-of- 
mass energy of y/s = 7TeV with the LHCb experiment at CERN. The detector (8) is 
a single- arm forward spectrometer covering the pseudorapidity range 2 to 5, designed 




r(g°(t) -» J/jjK") - T(B°(t) -» J/jjK° s ) 
T(B°(t) ->■ J/^KI) + T(B°(t) ->■ J/il)K%) 
Sjftiq sin(Am d t) - Cj^ K o cos(Am d t). 
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for the study of particles containing b or c quarks. It includes a high precision tracking 
system consisting of a silicon-strip vertex detector surrounding the pp interaction region, 
a large-area silicon-strip detector located upstream of a dipole magnet with a bending 
power of about 4 Tm, and three stations of silicon-strip detectors and straw drift-tubes 
placed downstream. The combined tracking system has a momentum resolution Ap/p that 
varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, and an impact parameter resolution 
of 20 \xm for tracks with high transverse momentum. Charged hadrons are identified 
using two ring-imaging Cherenkov detectors. Photon, electron and hadron candidates are 
identified by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic and a hadronic calorimeter. Muons are identified by a system composed 
of alternating layers of iron and multiwire proportional chambers. 

The analysis is performed on events with reconstructed B° — > J/tpK® candidates with 
subsequent J/ip — > and K® — > -n + n~ decays. Events are selected by the trigger 

consisting of hardware and software stages. The hardware stage accepts events if muon 
or hadron candidates with high transverse momentum (px) with respect to the beam 
axis are detected. In the software stage, events are required to contain two oppositely- 
charged particles, both compatible with a muon hypothesis, that form an invariant mass 
greater than 2.7GeV/c 2 . The resulting J/ip candidate has to be clearly separated from the 
production vertex (PV) with which it is associated on the basis of the impact parameter. 

Further selection criteria are applied offline to decrease the number of background 
candidates. The J/0 candidates are reconstructed from two oppositely-charged, well 
identified muons with p^ > 500 MeV/ c that form a common vertex with an invariant mass 
in the range 3035-3160 MeV/ c 2 . It is required that the J/0 candidate fulfils the trigger 
requirements described above. The K® candidates are formed from two oppositely-charged 
pions, both with (long K® candidate) or without (downstream K® candidate) hits in 
the vertex detector. Each pion must have p > 2 GeV/c and a clear separation from any 
PV. Furthermore, they must form a common vertex with an invariant mass within the 
range 485.6-509.6 MeV/ c 2 (long K° candidates) or 476.6-518.6 MeV/c 2 (downstream K° 
candidates). Different mass windows are chosen to account for different mass resolutions 
for long and downstream K® candidates. The K® candidate's decay vertex is required to 
be significantly displaced with respect to the associated PV. 

The B° candidates are constructed from combinations of J/0 and candidates that 
form a vertex with a reconstructed mass Trij^ K o in the range 5230-5330 MeV/ c 2 . The 
value of mj/^ K o is computed constraining the invariant masses of the /i + yU~ and 7r + 7r~ to 
the known J/0 and K® masses (9|, respectively. As most events involve more than one 
reconstructed PV, B° candidates are required to be associated to one PV only by requiring 
a clear separation from other PVs. Additionally, a significant decay time of the K® with 
respect to the B° decay vertex is required. 

The decay time t of the B° candidates is determined from a vertex fit to the whole 
decay chain under the constraint that the B° candidate originates from the associated 



PV 10 . Only candidates with a good quality vertex fit and with 0.3 < t < 18.3 ps are 
retained. In case more than one candidate is selected in an event, that with the best 
vertex fit quality is chosen. The fit uncertainty on t is used as an estimate of the decay 



2 



time resolution a t , which is required to be less than 0.2 ps. Finally, candidates are only 
retained if the flavour tagging algorithms provide a prediction for the production flavour 
of the candidate, as discussed in Section [3| 

Simulated samples are used for cross-checks and studies of decay time distributions. 
For the simulation, pp collisions are generated using Pythia 6.4 (TT] with a specific LHCb 



configuration 12 



Decays of hadronic particles are described by EvtGen 13 in which 

Til 



final state radiation is generated using Photos 
particles with the detector is implemented using the Geant4 toolkit 
Ref. [16 



The interaction of the generated 
as described in 
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3 Flavour tagging 

A mandatory step for the study of CP violating quantities is to tag the initial, i.e. 
production, flavour of the decaying B° meson. Since b quarks are predominantly produced 
in bb pairs in LHCb, the flavour tagging algorithms used in this analysis [17) reconstruct 
the flavour of the non-signal b hadron. The flavour of the non-signal b hadron is determined 
by identifying the charge of its decay products, such as that of an electron or a muon 
from a semileptonic b decay, a kaon from a b — > c — > s decay chain, or the charge of its 
inclusively reconstructed decay vertex. The algorithms use this information to provide a 
tag d that takes the value +1 (—1) in the case where the signal candidate is tagged as an 
initial B° (B°) meson. 

A careful study of the fraction of candidates that are wrongly tagged (mistag fraction) is 
necessary as the measured asymmetry is diluted due to the imperfect tagging performance. 
The mistag fraction (w) is extracted on an event-by-event basis from the combined per- 
event mistag probability prediction rj of the tagging algorithms. On average, the mistag 
fraction is found to depend linearly on r] and is parameterised as 

wfo) =Pi • (V~ (v)) +Po ■ (2) 

Using events from the self-tagging control channel B + — > J/ipK + , the parameters are 
determined to be p x = 1.035 ± 0.021 (stat) ± 0.012 (syst), p = 0.392 ± 0.002 (stat) ± 



0.009 (syst) and (rj) = 0.391 18 . The systematic uncertainties on the tagging calibration 
parameters are estimated by comparing the tagging performance obtained in different decay 
channels such as B° — > J/i/jK*°, in B + and B~ subsamples separately, and in different 
data taking periods. 

The difference in tagging response between B° and B° is parameterised by using 

w = uj(t}) ± , (3) 

where the + (— ) is used for a B° (B°) meson at production and Ap Q is the mistag fraction 
asymmetry parameter, which is the difference of po for B° and B° mesons. It is measured 
as Apo = 0.011 ± 0.003 using events from the control channel B + — > J/ipK + . By using 
Apo in the analysis, the systematic uncertainty on the po parameter is reduced to 0.008. 
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The effect of imperfect tagging is the reduction of the statistical power by a factor e tag V 2 , 
where e tSLg is the tagging efficiency and V = 1 — 2w is the dilution factor. The effective e tag 
and V values are measured as e ta , g = (32.65 ± 0.31)% and V = 0.270 ± 0.015, resulting 
in Eta.gV 2 = (2.38 ± 0.27)%, where combined systematic and statistical uncertainties are 
quoted. 

4 Decay time acceptance and resolution 

The bias on the decay time distribution due to the trigger is estimated by comparing 
candidates selected using different trigger requirements. In the selection, the reconstructed 
decay times of the B° — > J/ipK® candidates are required to be greater than 0.3 ps. This 
requirement makes the acceptance effects of the trigger nearly negligible. However, some 
small efficiency loss remains for small decay times. Neglecting this efficiency loss is treated 
as a source of systematic uncertainty. 

A decrease of efficiency is also observed at large decay times, mostly affecting the 
candidates in the long K® subsample. This can be described with a linear efficiency 
function with parameters determined from simulated data for the downstream and long 
Kg subsamples separately. The efficiency function is then used to correct the description 
of the decay time distribution. 

The finite decay time resolution of the detector leads to an additional dilution of the 
experimentally accessible asymmetry. It is modelled event-by-event with a triple Gaussian 
function, 

nt - ty t ) = e u^- ex P ( Jt ~J~ b ° t)2 ) , (4) 

where t is the reconstructed decay time, t' is the true decay time, and Oj is the per-event 
decay time resolution estimate. The parameters are: the three fractions fi, which sum to 
unity, the three scale factors s i: and a relative bias b, which is found to be small. They 
are determined from a fit to the t and a t distributions of prompt J/0 events that pass the 
selection and trigger criteria for J?°— >• J/ipK®, except for decay time biasing requirements. 
The parameters are determined separately for the subsamples formed from downstream 
and long K® candidates. This results in an average effective decay time resolution of 55.6 fs 
(65.6 fs) for candidates with long (downstream) Kg. 

5 Measurement of Sj^ K o and Cj^ K o 

The analysis is performed using the following set of observables: the reconstructed mass 
rrij^ K o, the decay time t, the estimated decay time resolution a t , the flavour tag d, and 
the per-event mistag probability 77. The CP observables Sj^ K o and Cj/^ K o are determined 
as parameters in an unbinned extended maximum likelihood fit to the data. 
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Due to different resolution and acceptance effects for the downstream and long K® 
subsamples, a simultaneous fit to both subsamples is performed. In each subsample, the 
probability density function (PDF) is defined as the sum of two individual PDFs, one for 
each of the components of the fit: the B° signal and the background. The latter component 
contains both combinatorial background and mis-reconstructed 6-hadron decays. 

The reconstructed mass distribution of the signal is described by the sum of two 
Gaussian PDFs with common mean but different widths. Only the mean is shared between 
the two subsamples. The background component is parameterised as an exponential 
function, different for each subsample. 

The signal and background distributions of the per-event mistag probability 77 are 
modelled with PDFs formed from histograms obtained with the sPlot technique [19] on 
the reconstructed mass distribution. In both subsamples the same signal and background 
models are used. 

The distributions of the estimated decay time resolution at are different in each 
component and each subsample. Hence, no parameters are shared between subsamples or 
components. All a t PDFs are modelled with lognormal functions 



where M at is the median and k the tail parameter. The background components in both 
subsamples are parameterised by single lognormal functions. For the signal a sum of 
two lognormals with common (different) median parameter (s) is chosen for the long K® 
(downstream K®) subsample. 

The background PDFs of the decay time are modelled in each subsample by the sum of 
two exponential functions. These are convolved with the corresponding resolution function 
lZ(t — t'\at). The parameters are not shared between the two subsamples. The background 
distribution of tags d is described as a uniform distribution. 

The signal PDF for the decay time simultaneously describes the distribution of tags d, 
and is given by 

V(t,d\a t ,v) =e{t)-Vcp{t',d\a t ,Tj)®n{t-t'\at) , (6) 

with 

V C p(t', d\a u rj) oc e" 1 ' /r (l - dA Po - dA P (l - 2u{q)) 

- (d(l - 2u)(r))) - A P (1 - dAp ))S J/i}K o sinAm d t' 

+ - 2u>(rj)) - A P (1 - dAp Q ))C J/i>K o cos Am/) . (7) 

This PDF description exploits time-dependent asymmetries, while its normalisation adds 
sensitivity by accessing time-integrated asymmetries. The lifetime r, the mass difference 
Airid, and the CP parameters Sj/^ K o and Cj^ K o are shared in the PDFs of the downstream 
and long K® subsamples, as well as the asymmetry Ap = (R^o — Rb°)/(Rb° + Rb°) of 
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the production rates R for B° and B° mesons in pp collisions at LHCb. The latter value 
has been measured in Refs. 



20,21 



to be A P = -0.015 ±0.013. 
In the fit all parameters related to decay time resolution and acceptance are fixed. The 
tagging parameters and the production asymmetry parameter are constrained within their 
statistical uncertainties by Gaussian constraints in the likelihood. The fit yields 



S 



0.73 ±0.07 



0.03 ±0.09, 



with a correlation coefficient p{Sj^ko,Cj^j(o) = 0.42. Both of the uncertainties and the 
correlation are statistical only. The lifetime is fitted as r = 1.496 ± 0.018 ps and the 
oscillation frequency as Am d = 0.53 ± 0.05 ps -1 , both in good agreement with the world 
averages (7 



22 



The mass and decay time distributions are shown in Fig. [I] The measured 



signal asymmetry and the projection of the signal PDF are shown in Fig. [2] 




Figure 1: Invariant mass (left) and decay time (right) distributions of the B° — > J/if)K® 
candidates. The solid line shows the projection of the full PDF and the shaded area the 
projection of the background component. 



6 Systematic uncertainties 

Most systematic uncertainties are estimated by generating a large number of pseudo- 
experiments from a modified PDF and fitting each sample with the nominal PDF. The 
PDF used in the generation is chosen according to the source of systematic uncertainty 
that is being investigated. The variation of the fitted values of the CP parameters is used 
to estimate systematic effects on the measurement. 

The largest systematic uncertainty arises from the limited knowledge of the accuracy of 
the tagging calibration. It is estimated by varying the calibration parameters within their 
systematic uncertainties in the pseudo-experiments. Another minor systematic uncertainty 
related to tagging emerges from the difference of tagging efficiencies of B° and B° which 
is measured as Ae tag = 0.000 ± 0.001 in B + J/ipK + . 
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Figure 2: Time-dependent asymmetry (Ng — N B o)/(N^ + N B o). Here, N B o (iVg ) is 
the number of 5° — >■ J/ipK® decays with a Z?° (B°) flavour tag. The data points are 
obtained with the sPlot technique, assigning signal weights to the events based on a fit to 
the reconstructed mass distributions. The solid curve is the signal projection of the PDF. 
The green shaded band corresponds to the one standard deviation statistical error. 



The effect of an incorrect description of the decay time resolution model is derived from 
pseudo-experiments in which the scale factors of the resolution model are multiplied by a 
factor of either 0.5 or 2 in the generation. As the mean decay time resolution of LHCb is 
much smaller than the oscillation period of the B° system this variation leads only to a 
small systematic uncertainty. The omission of acceptance effects for low decay times is 
estimated from pseudo-experiments where the time-dependent efficiencies measured from 
data are used in the generation but omitted in the fits. Additionally, a possible inaccuracy 
in the description of the efficiency decrease at large decay times is checked by varying the 
parameters within their errors, but is found to be negligible. 

The uncertainty induced by the limited knowledge of the background distributions is 
evaluated from a fit method based on the sPlot technique. A fit with the PDFs for the 
reconstructed mass is performed to extract signal weights for the distributions in the other 
observable dimensions. These weights are then used to perform a fit with the PDF of 
the signal component only. The difference in fit results is treated as an estimate of the 
systematic uncertainty. 

To estimate the possible bias in the CP parameters emerging from the fit method itself, 
the method is probed with a large set of pseudo-experiments. The bias is found to be 
0.004 for Sj/^ko and 0.005 for Cj/^ K o. These values are added as systematic uncertainties. 
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Table 1: Summary of systematic uncertainties on the CP parameters. 



Origin 






Tagging calibration 


0.034 


0.001 


Tagging efficiency difference 


0.002 


0.002 


Decay time resolution 


O.OOf 


0.002 


Decay time acceptance 


0.002 


0.006 


Background model 


0.012 


0.009 


Fit bias 


0.004 


0.005 


Total 


0.036 


0.012 



The uncertainty on the scale of the longitudinal axis and on the scale of the momen- 
tum 23 sum to a total uncertainty of < 0.1% on the decay time. This has a negligible 



effect on the CP parameters. Likewise, potential biases from a non-random choice of the 
B° candidate in events with multiple candidates are found to be negligible. 

The sources of systematic effects and the resulting systematic uncertainties on the CP 
parameters are quoted in Table [T] where the total systematic uncertainty is calculated by 
summing the individual uncertainties in quadrature. 

The analysis strategy makes use of the time-integrated and time-dependent decay rates 
of B° — > J/ip K® decays that are tagged as B°/B° meson. Cross-check analyses exploiting 
only the time-integrated or only the time-dependent information show that both give 
results that are in good agreement and contribute to the full analysis with comparable 
statistical power. 



7 Conclusion 

In a dataset of l.Ofb -1 collected with the LHCb detector, approximately 8200 flavour 
tagged decays of B° — > J/ipK® are selected to measure the CP observables Sjj^ K o and 
Cj/^ko, which are related to the CKM angle (3. A fit to the time-dependent decay rates of 
B° and B° decays yields 

S J/i}K o = 0.73 ± 0.07 (stat) ± 0.04 (syst), 
C J/i>K o = 0.03 ± 0.09 (stat) ± 0.01 (syst), 

with a statistical correlation coefficient of p(Sj/^ K o,Cj/^ K o) = 0.42. This is the first 



significant measurement of CP violation in B° — > J/ip K® decays at a hadron collider 24 
The measured values are in agreement with previous measurements performed at the B 
factories [5j|6] and with the world averages [7j. 
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